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INTRODUCTION 


It  is  typically  of  great  Interest  and  practical  importance  to  accelerate 
a  vibration  or  thermal  cycling  qualification  or  acceptance  test  in  the 
laboratory  from  the  actual  service  conditions.  The  test  duration  is 
compressed  by  a  relatively  large  factor  (e.g.  1000)  with  an  attendant 
increase  in  applied  vibration  or  temperature  limit  levels.  The  relationship 


between  test  duration  compression  and  level  Increase  th 
by  the  industry  is  as  follows : 


\has  been  used 


Cl) 


where  T^  -  service  life 
T£  »  test  duration 
-  service  level 
Lj  *  test  level 
a  -  constant 


Various  industry  groups  assign  various  values  to  o  for  a  variety  of  reasons 
and  conservativeness.  Values  of  a  between  4  to  9  have  been  used.  There 
is  much  disagreement  between  groups  on  the  value  assigned  to  a.  As  an 
example  for  a  •  4  and  a  1000  hour  service  life  compressed  into  a  one  hour 
test 

-2-  -  (1000) 1/4  -  5.6 

L1 

Thus,  the  "laboratory"  test  level  should  be  5.6  larger  than  the  service 
level  in  order  that  the  cumulative  fatigue  damage  is  the  same  in  both  cases. 
Fatigue  damage  is  not  to  be  misinterpreted  as  fatigue  failure  (i.e.  fracture). 
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Fatigue  failures  may  or  may  not  occur  depending  upon  the  magnitude  of 
the  cumulated  damage.  The  previous  equation  can  be  interpreted  to  mean 
that  the  potential  for  fracture  will  be  the  same  for  a  service  duration 
Tj  at  level  as  for  a  test  duration  T2  at  Level  L2  . 

This  paper  shows  a  logical  basis  for  selecting  the  form  of  the  relationship 
and  assigning  values  to  the  parameters.  For  Fatigue  (i.e.  no  initial  flaws) 
the  form  of  the  correct  expression  is  as  previously  shown.  The  introduction 
of  Fracture  Mechanics  effects  (i.e.  initial  flaws)  results  in  a  different 
form  of  the  expression  that  relates  levels  ,  durations  and  different 
parameter  values.  Non-linear  damping  effects  are  also  included.  Conditions 
of  similitude  between  the  service  and  accelerated  test  environments  are 
identified.  Such  conditions  must  be  ensured  if  the  developed  equations  are 
to  be  applied  appropriately  and  accurately  (e.g.  no  new  failure  mechanism 
should  be  introduced  at  Che  accelerated  test  level). 

Techniques  are  described  that  allow  compensating  during  the  accelerated 
test  for  differences  in  response  stress  spectra  or  distribution  of  stress 
peaks  (clipping)  between  the  service  and  laboratory  environments. 
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APPROACH  SUMMARY 


*>1 


(3) 


Dx  “  damage  cumulated  in  service  environment 

Nx  *  number  of  applied  stress  cycles  at  stress  level  AS^  , 

Nf^  ■  number  of  stress  cycles  to  failure  at  stress  level  1 
Do  -  damage  in  accelerated  environment 

Nx  *  number  of  applied  stress  cycles  at  stress  level  AS2  1  o ^ 
Nf^  “  number  of  stress  cycles  to  failure  at  stress  level  2 
AS  *  sinusoidal  stress  range,  peak-peak  stress 
a  •  random  stress  rms  level 

Fatigue  failure  (fracture)  occurs  when  or  I>2  *  1. 


N1  *  flTl  »  n2  “  f2T2 

f2  -  frequency  of  stress  cycles  at  level  AS^  ,  o^ 

«  rim*  or  duration  of  applied  stress  at  level  AS^ 
f2  ■  frequency  of  stress  cycles  at  level  ASj  ,  Ox 
Tj  “  time  or  duration  of  applied  stress  at  level  ASj 


Time  compression  ratio 


for  f^ 


al 

a2 


(4) 


n2  <  nl 

Accelerated  test  level  ■  £12  ,  12.  (5) 

ASx 

Fatigue  failures  may  or  may  not  occur  at  either  the  environmental  or 
accelerated  test  levels.  Whether  they  occur  is  not  of  interest  for  this 
analysis.  It  is  of  interest,  however,  that  the  accumulated 
fatigue  damage  be  the  same  for  both  conditions  (i.e.  ■  D2).  For  example. 


if  602  of  life  (no  fatigue  failure)  is  accumulated  after  cycles  at  the 
service  environmental  stress  level  AS^  or  it  is  desired  to  find  the 
accelerated  test  stress  level  AS 2  or  O2  that  will  correspondingly  accumulate 
the  same  60%  of  life  after  only  N2  cycles  (i.e.  •  D2  ■  0.60).  Also  if 

120%  of  life  (fatigue  failure)  is  accumulated  at  the  service  level,  it  is 
desired  to  accumulate  120%  of  life  at  the  accelerated  test  level  (i.e.  ■ 

D2  *  1.20). 

Thus 


“1 

Ni_ 

Nfi 


u2 
N- 

4 


(6) 

(7) 


or 


(8) 


The  fatigue  and  test  parameters  are  related  directly  to  and  , 
not  to  and  N2.  The  analysis  that  follows,  therefore,  will  be  in 
terms  of  N^/Nf^.  This  ratio  is  the  same  as  ^/N^.  Expressions  in  terms 
of  fatigue  failure  parameters  are  a  mathematical  and  engineering 
necessity  but  should  not  be  misinterpreted  to  mean  that  fatigue  failures 
will  occur  at  either  level. 


The  ratio  (i.e.  the  ratio  ^/N^)  will  be  a  function  of  the 

corresponding  stress  levels.  This  equation  will  then  be  rearranged  so 
that  the  ratio  AS2/AS2  or  can  be  solved  for  in  terms  of  the  "time" 

compression  factor  N^/^. 


The  fatigue  process  can  be  characterized  by  crack  initiation,  stable  crack 
propagation  and  fracture  ( i.e .  unstable  crack  growth  when  the  crack  size 
equals  the  critical  crack  si2e  value).  Each  of  these  individual  processes 
is  directly  a  function  of  the  stress  or  strain  level  and  the  number 
of  applied  stress  cycles.  The  fatigue  process  is  only,  indirectly  related 
to  the  input  vibration  acceleration  level  to  the  structural  elements  in, 
say,  an  electronic  "black  box"  or  to  the  mission  or  test  duration. 

Consequently  relationships  have  been  established  that  relate  stress  levels 
to  input  vibration  acceleration  levels  and  number  of  stress  cycles  to 
mission  or  test  duration.  This  allows  the  input  vibration  acceleration 
levels  x_/x.  to  be  functionally  related  to  the  time  compression  factor 

Ti/t2. 

Linear  and  non-linear  dependence  of  stress  upon  input  vibration  level  is 
included  in  the  term  r>.  n  ■  1  corresponds  to  a  linear  relationship.  It 
is  shown  that  n  -  0.714  for  sinusoidal  vibration  inputs  and  n  »  0.833  for 
random  vibration  inputs  where  the  predominant  damping  mechanism  is  internal 
stress-strain  hysteresis  damping.  For  n  <  1  the  effective  damping  at 
resonance  increases  more  than  proportionally  with  an  increase  in  input 
vibration  level,  n  >  1  applies  to  cases  where  the  predominant  damping  mechanism 
is  Coulomb  friction  or  where  the  effective  spring  stiffness  increases  with 
input  vibration  level,  as  examples. 

Fracture  Mechanics  effects  (i.e.  initial  cracks  or  flaws)  have  been  distinguished 
from  the  usual  Fatigue  effects  (i.e.  no  initial  flaws).  The  major  difference 
is  in  the  number  of  stress  cycles  (hence,  time)  required  to  initiate  cracks 
for  Fatigue.  C racks  (either  actual  or  postulated)  already  exist  for  Fracture 
Mechanics. 
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The  solution  of  the  accelerated  stress  level  ratio  in  terms  of  the  stress 
cycle  "compression  ratio"  and  the  accelerated  vibration  input  level  ratio 
in  terms  of  the  time  compression  ratio  involves  solving  transcendental 
functions.  Computer  program  in  Basic  Language  that  accomplish  that  task 
are  included. 


CONDITIONS  OF  SIMILITUDE 


Certain  conditions  of  similitude  must  be  imposed  upon  the  service  and 
laboratory  accelerated  test  environments  if  the  developed  mathematical 
relationships  are  to  be  appropriately  and  accurately  applied.  The 
fundamental  hypothesis  is  that  the  damage  states  and  damage  rates  must 
be  the  same  for  both  environments.  Specifically  the  states  of  stress 
(torsion,  bending,  axial),  the  corresponding  strengths,  the  resonant 
mode  shapes,  the  internal  response  stress  spectrum  shapes,  the  stress 
peak  distribution,  and  the  type  and  location  of  failure  mechanisms  must 
be  the  same  for  both  environments. 

Differences  in  temperature,  rate  of  stressing,  corrosive  environments,  and 
other  environmental  effects  (e.g.  "purple  plague"  that  can  result  from 
combined  high  temperature  and  humidity)  between  the  service  and  test 
conditions  may  violate  conditions  of  similitude  for  some  structural  elements. 
Violation  will  occur  if  the  above  factors  are  sufficient  to  alter  the 
material's  fatigue  strength  (i.e.  fatigue  curve  parameters)  between  the  two 
environments. 

Threshold  sensitive  or  ocher  non-linear  response  effects  in  general  tend 
to  violate  conditions  of  similitude.  In  some  cases  lack  of  similitude  can 
be  quantitatively  compensated  for.  Severe)  examples  are  included  in 
Appendix  A. 

The  condition  that  the  shape  of  the  vibration  acceleration  input  spectra 
must  be  the  same  for  both  environments  has  purposely  been  omitted  from  the 
previously  listed  conditions.  This  is  because  the  fatigue  damage  state 
and  race  are  only  indirectly  related  to  the  input  acceleration  spectrum. 
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They  are  directly  related  to  the  response  stress  spectrum  at  the  location 
where  damage  is  accumulating.  Stress  Is  herein  defined  as  the  internal 
force  per  unit  area  in  a  material  that  results  from  the  application  of  an 
external  load.  The  input  vibration  acceleration  to  a  "black  box"  is 
defined  as  the  kinematic  motion  response  at  the  load  transfer  path  input 
location  that  results  from  applied  vibratory  loads  to  the  black  box  and 
adjacent  structural  members.  The  input  vibration  acceleration  is  not  an 
"applied  stress”  using  the  above  definitions. 


S 


FATIGUE 


SINE  VIBRATION 


From  reference  [1]  a  material's  sine  fatigue  curve  Is 


n  *1  for  linear  damping  (refer  to  Appendix  B) 

x2  ■  represents  the  accelerated  test  input  acceleration  level  (g's) 
x^  -  service  vibe  Input  acceleration  level  (g's) 


Tp  -  corresponding  durations 


FATIGUE 


RANDOM  VIBRATION 


From  reference  [ 1]  the  random  fatigue  curve  equation  is 
o  »  C  N'1/6  (13) 

where 

a  •  rms  stress  (KSI) 

C  »  constant  (KSI) 

N  *  average  cycles  to  failure 

B  ■  slope  parameter  of  the  material's  sine  fatigue  curve 


N  *  fT 


(1A) 


where 

f  -  center  frequency  of  narrow-band  response  (Hz) 

T  -  duration 


In  a  fashion  similar  for  the  sinusoidal  case  the  accelerated  test  level  is 


3  -ftp  ft) 


1/8 


(15) 


CA  **n 


(16) 
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<J  (17) 

n  *  1  for  linear  damping 

•  4 

■  service  vibration  input  acceleration  rms  level 

»  duration  of  service  vibe 

•• 

»2  •  accelerated  test  vibe  input  acceleration  rms  level 

T2  “  duration  of  accelerated  test 


The  above  equations  assume  that  the  ratio  of  the  input  Power  Spectral  Density 
(PSD)  W0  in  the  vicinity  of  resonance  is 


(18) 


Otherwise,  the  conditions  of  similitude  will  be  violated. 
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LOW  CYCLE  FATIGUE 


Mechanical  loads  or  deformations  that  are  of  sufficient  magnitude  to 
stress  the  material  into  the  plastic  (i.e.  inelastic)  region  of  its 
stress-strain  curve  are  associated  with  short  fatigue  lifes.  This  is  often 
referred  to  as  low  cycle  fatigue.  Fatigue  lifes  typically  extend  up  to 

A 

approximately  10  stress  cycles.  The  high  cycle  fatigue  (i.e.  elastic 

4 

stress-strain)  typically  extends  beyond  10  cycles. 

The  Coffin-Manson  low  cycle  fatigue  expression  relates  the  applied  plastic 
strain  amplitude,  the  material's  ductility  and  cycles  to  failure  for 
cyclically  induced  strains  in  mechanical  systems.  Static  stresses  do  not 
affect  fatigue  life  in  the  low  cycle  region  and  therefore  can  be  ignored 
for  this  analysis. 


H-  ■  «'f  «»£> 


■  applied  plastic  strain  amplitude  (in/in) 


■  fatigue  ductility  coefficient  (in/in) 


Nf  **  cycles  to  failure 


6  “  slope  parameter  *  2  for  most  structural  materials 


The  above  strains  are  "true"  strains  which  include  changes  in  the 
strained  specimen's  cross-sectional  area  under  load  as  compared  to 
"engineering"  strains  which  are  based  upon  the  specimen's  elongation 
relative  to  its  original  length.  Until  specimen  necking  occurs: 


** 


'  true 


in  (1  + 


eeng'  g'1 


(20) 


"Engineering  strain  is  usually  more  convenient  to  use  than  "true"  strain. 
The  Coffin-Manson  can  be  modified  I 2J  to  give  > 

-1/6 

e.  -  eu  (2  Hf)  (21) 

e  ■  applied  "engineering"  plastic  strain  amplitude  (in/ln) 
eu  *  material  ductility;  ultimate  percent  elongation  (in/in) 


Nj  and  6  are  the  same  as  before 


TABLE  I  CORRESPONDING  STRAIN  PARAMETERS 


Eng'g 

True 

e 

Ae/2 

i 

e 

u 

e  f 

Using  true  strains  the  service  and  accelerated  test  levels  can  be  related 
as  follows : 


(22) 


where  the  subscript  1  applies  to  the  service  environment  and  the 
subscript  2  applies  to  the  accelerated  test  environment. 


Engineering  strains  can  be  substituted  into  the  above  expression.  The 
results  will  then  be  accurate  if  no  necking  of  the  structural  element  occurs. 
The  results  will  be  conservative  if  necking  does  occur,  because  the  material 
is  actually  more  ductile  than  given  credit  for. 
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The  above  expression  applies  to  all  forms  of  cyclic  strain.  Strains 
resulting  from  temperature  cycling  is  typical. 


As  an  example  consider  a  glass  epoxy  multi-layer  board  (MLB)  that  is  to 
be  subjected  to  temperature  cycling.  The  differential  expansion  rate  in 
a  direction  perpendicular  to  the  plane  of  the  board  between  the  epoxy 
and  the  electrodeposited  copper  plated-through-holes  (PTH)  is  non-linearly 
related  to  temperature.  Assume  that  the  service  temperature  cycle  limits 
of  0°C  to  +95°C  produces  an  applied  strain  amplitude  of  8.45  x  10~4  in/in 
in  the  middle  region  of  the  PTH's  where  there  is  a  potential  for  circumferential 
cracking.  Assume  that  the  quantity  of  service  temperatures  cycles  is  7000 
cycles.  It  is  desired  to  find  the  accelerated  test  temperature  range  to 
cumulate  the  same  fatigue  damage  in  only  2550  cycles. 


c 2  *  0.0014  in/in 

Measured  strain  amplitude  versus  temperature  limits  indicates  that  0.0014 
in/in  corresponds  to  limits  of  -65 °C  to  +125 "C  for  the  particular  MLB. 

Care  must  be  taken  to  ensure  that  a  new  failure  mechanism  is  not  introduced 


(e.g.  PTH  comer  cracks). 


Care  muse  be  exercised  when  considering  materials  whose  fatigue  properties 
are  rate  or  test  temperature  sensitive  £3iJ.  Consider  63  -  37  Tin-Lead  Solder 
plastically  stressed  in  reversed  shear; 


TABLE  II  SOLDER  SHEAR  FATIGUE  PARAMETERS 


TEST 

TEMPERATURE 

(C°) 

SHEAR 

STRAIN 

RATE 

(cycles  per 
minute) 

e 

FATIGUE 
EXPRESSION 
(in /in) 

25 

1/15 

Ac  *  0.531  n~°* 381 

25 

5 

Ac  -  0.560  tT0’302 

100 

5 

Ac  -  0.488  n“*0348 

15 
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fracture  mechanics  effects 


The  primary  fracture  mechanic's  effects  are  those  due  to  initial  cracks 
(flaws)  that  are  either  actual  or  hypothesized.  Such  cracks  reduce  fatigue 
life.  They  either  exist  in  the  structural  material  as  metallurgical 
inclusion  or  dislocations  or  are  introduced  during  manufacturing  fabrication 
and  assembly  operations.  They  can  also  be  created  by  temporary  overloads 
into  the  plastic  stress  regions. 


FATIGUE  CURVES 


From  typical  fatigue  curve  data  (e.g.  reference  £4?  or  reference  [5})where 
any  initial  flaw  sizes  are  approximately  zero  the  usual  form  of  the  fatigue 
curve  is 


-  A  IT1^  . 


For  7075-T6  Aluminum  Alloy 


AS 

2 


-0.104 


180  N  KSI 

From  Fracture  Mechanics  (See  Appendix  C) 
1.698  x  10 


Nf  - 


AS 


l  fjL  - 

_ai 


-3  2 

7.83  x  10  AS 


2 


The  above  Nf  expression  applies  for  a  particular  geometry. 


Both  types  of  fatigue  curves  are  plotted  in  figure  1  .  It  can  be  seen 

that  even  a  small  value  of  aj^  reduces  fatigue  life.  Further  observations 

are: 

a)  The  slope  parameter  of  the  usual  fatigue  curve  (  a^  *  0  )  is 
8  -  9.65  for  7075-T6. 
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b)  The  slope  parameter  of  the  fracture  mechanics  fatigue  curve 


(  ai  -  0  )  is  0.  6  -  4  for  7075-T6. 

NOTE:  6  f  8 

c)  For  large  Nf  (7075-T6) 


I 


ACCELERATED  SINE  STRESS  LEVEL 


Appendix  D'  (FRACTURE  MECHANICS  ACCELERATED  SINUSOIDAL  TEST  STRESS  LEVEL) 
shows  the  derivation  of  the  accelerated  sine  test  stress  level  AS2 
given  the  environmental  stress  level  AS^  and  the  corresponding  test 
duration  (i.e.  cycles)  N^  and  N^  .  The  transcendental  function  cannot 
be  normalized  (e.g.  AS2/AS2  versus  ^/N^)  because  of  the  inherent 
non-linearities  in  the  fracture  mechanics  correction  factor  X.  Thus, 

AS1  and  must  be  assigned  specific  values. 

The  Basic  Language  computer  program  PL-2  solves  for  AS2 ,  the  accelerated 
stress  level.  The  inputs  are  N^  ,  N2  ,  8  ,  Y  ,  a^.  ,  AKC  ,  AS^  .  The 
listing  shows  typical  parameter  values  for  7075-T6  Aluminum  Alloy.  It 
should  be  noted  that  AKC  was  chosen  to  be  20  KSI^IN.'  £5]  This  value  was 
the  lowest  (hence,  the  most  conservative)  value  published  in  the 
literature.  Much  higher  values;  unfortunately,  have  also  been  published. 
Therefore,  care  must  be  exercised  in  using  published  data. 

It  is  recommended  that  the  most  reliable  data  is  that  obtained  using  test 

method  ASTM  E647-78T,  "Tentative  Test  Method  for  Constant  Load  Amplitude 

—8 

Fatigue  Crack  Growth  Rates  Above  10  m/ cycle."  This  method  uses  as  an 

applicable  document  ASTM  399  "Test  for  Plane-Strain  Fracture  Toughness  of 
Metallic  Materials."  Reference  C 6 ]  uses  this  test  method. 

Several  runs  using  PL-2  are  included  for  7075-T6.  See  Groups  I  -  III. 
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Figure  2  shows  Che  converted  results.  The  accelerated  levels  are  not 
sensitive  to  a^  values  for  small  time  compression  ratios  but  are 
sensitive  for  large  time  compression  ratios. 
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FIGURE  Z  ACCELERATED  5INE  VI&E  TEST 

LEVEL  FOR  0.£  =  0.007. 0.05*  .O.f  INCHES 


ACCELERATED  SINUSOIDAL  VIBE  TEST  ACCELERATION  LEVEL 


In  che  previous  section  Che  accelerated  stress  &S2  was  calculated 
given  AS^  ,  N-j_  and  N2  •  The  corresponding  sine  acceleration  -  time 
parameters  can  be  calculated  from  the  stress  -  cycles  parameters.  Use 
PL-2  to  obtain  the  stress-cycle  parameters.  Then  compute  the  acceleration¬ 
time  parameters  using  the  following  equations. 


(26) 


(27) 


T1  “  Vfn  (28) 

T2  m  Vfn  (29) 


where  fQ  ■  resonant  frequency 


2  ' 


See  Appendix  E  (SINE-RANDOM  EQUIVALENCE)  for  an  example  of  C 


ACCELERATED  RANDOM  VIBRATION  TEST  LEVEL 

Appendix  F  (Derivation  0f  Accelerated  Random  Vibe  Test  Acceleration  Level) 

shows  the  derivation  of  the  accelerated  random  vibe  test  acceleration  input 

*• 

level  *2  (g  rms)  given  the  environmental  input  acceleration  level  x^ 

(g  rms)  and  the  corresponding  test  durations  T^  and  T2  .  The  expression  is: 


This  transcendental  function  cannot  be  normalized  in  closed  form  in  terms  of 

M  M 

x2'xl  versus  T1/T2  because  of  the  inherent  non-linearities  in  the  fracture 
mechanics  correction  factor  X  .  Thus  x^  and  T^  must  be  assigned 
specific  values. 

The  Basic  Language  Computer  program  listing  PL-3  solves  the  previous 
transcendental  function.  The  program  inputs  are  the  test  durations  T1  and 
T2,  6  ,  n  ,  ^  ,  C4  ,  A  ,  C  ,  &KC  and  Y  .  See  PL-3  for 

further  details.  The  program  output  is  the  accelerated  random  vibe  test  input 
rms  acceleration  value,  X2  . 

Figure  3  shows  a  plot  of  ”2  versus  T2  for  several  values  of  a^  . 

•  • 

Xj_  ■  1  g  rms  ,  T1  -  1000  hours  and  n  »  0.833  .  These  curves  are  almost 
straight  lines  on  a  log-log  plot.  Increasing  values  of  a^  from  0.007  inch 
to  0.100  inch  does  not  greatly  alter  the  curves  for  the  T2  values  shown. 

For  lower  values  of  T2  the  curves  deviate  more  from  each  other  similar  to 
those  of  figure  2 
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Figure  4  has  the  same  parameter  values  as  figure  3  except  n  ■  1  . 

The  results  are  similar.  Figure  5  has  n  ■  1.2  .  Again  the  results  are 
similar.  The  curves  of  all  three  figures  have  slopes  on  the  log-log  plots 
as  follows : 


slope 


_1_ 

n6 


Thus,  it  might  be  expected  that  the  results  should  be  sensitive  to  ri  values. 
Figure  6  confirms  this  expectation.  Since  6  is  a  material  property,  the 
results  are  also  sensitive  to  the  material.  Table  V  shows  several  material 
6  values . 


TABLE  V  6  Values 


MATERIAL 

e 

A-286 

3.24 

A  471  a  4 

1.4 

Cr-Mo-V 

4.09 

4340 

4.65 

7075-T6 

4.00 

It  will  be  noted  that  for  all  fracture  mechanic's  examples  in  this  study 
figure  7  will  apply  for  simplicity. 


tz  (  hours') 

i 

Figure  3  Accelerated  random  input  (>>  =  0.$ 33} 
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figture  4  Accelerated  random  input  = 


V*«. i , 


FIG-ORE  5  ACCELERATED  RANDOM  IKIPOT  (">7=1.2) 


FlGrURE  (,  ACCELERATED  RANpoM  INPUT  C  7  VARIABLE^ 


CENTER  CRACKED  STRIP 
LOADED  IM  TENSlOrt 


PL-3 


ACCELERATED  RANDOM  TEST 
LEVEL  COMPUTATION  PROGRAM 
(BASIC  LANGUAGE) 


LOW  CYCLE  FATIGUE 


The  previous  development  of  fracture  mechanics  effects  was  restricted  to 
elastic  stress  fields  where  the  use  of  the  stress  intensity  factor  AK 
and  fracture  toughness  AKC  are  well  established  and  readily  applied. 
The  crack  growth  rate  expression  used  was  that  of  Paris  17]: 


da 

dN 


co  AK® 


(31) 


where 


AK 


Y  ASa 


1/2 


(32) 


The  high  strain  fatigue  (i.e.  low  cycle  fatigue  region)  crack  growth  rate 
characterization  is  not  as  well  established. 

The  most  accurate  characterization  is  the  J-Integral  explored  by  Dowling 
and  Begley  [8]  [9]  [lO] .  J  is  a  line  integral. 


The  Dowling  and  Begley  expression  is: 


da 

dN 


C1  ajy 


(35) 


where  AJ  is  the  range  of  the  energy  line  integral  J  ,  and  C2  and  y 
are  material  dependent  constants.  This  expression  has  general  use  to  all 
materials  [ 12] . 


Mowbray  £ic0  has  shown  that  this  relationship  also  reduces  to  the  Coffin- 
Manaon  low  cycle  fatigue  expression.  An  important  aspect  of  the  Dowling 
and  Begley  work  is  that  only  the  loading  during  crack  face  opening  results 
in  damage. 
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At  present  there  is  one  objection  to  applying  the  J-Integral  to  fatigue 
crack  growth  and  that  pertains  to  the  mathematical  definition  of  J  Cll] . 

It  is  mathematically  valid  within  the  limits  of  deformation  plasticity 
theory,  which  precludes  unloading.  Dowling  approaches  this  objection  on 
the  basis  that  J  may  have  more  applicability  than  the  current  mathematical 
definition  indicates.  More  test  data  will  help  resolve  this  issue. 

There  is  difficulty  at  present  with  the  practical  application  of  the 
approach,  the  determinition  of  J  versus  crack  length  "a"  relationships  £ll3 . 
There  are  only  a  limited  number  of  configurations  for  which  J  is  known 
or  can  be  directly  measured.  However,  any  approach  involving  non-linear 
material  behavior  will  have  similar  difficulties. 

The  previously  developed  elastic  accelerated  test  level  equations  are  in 
terms  of  the  Paris  equation  parameters;  namely,  cQ  ,  AK  and  6  . 

Those  same  equations  can  be  used  to  determine  the  inelastic  accelerated 
test  levels  by  substituting  as  follows: 


TABLE  VI  EQUIVALENT  PARAMETERS 


IN-PLACE 

SUBSTITUTE 

OF 

DOWLING-BEGLEY 

AK 

AJ 

AK 

c 

AJC 

co 

C1 

e 

y 
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SINE-RANDOM  EQUIVALENCY 


Appendix  E  (SINE-RANDOM  EQUIVALENCE  DERIVATION)  derives  the  relationship 
between  the  sinusoidal  "black  box"  vibe  input  level  x  and  the  random 
vibe  input  power  spectral  density  Wq  in  the  vicinity  of  structural 
resonance  that  will  cumulate  the  same  fatigue  damage  in  the  same  test 
time.  The  desired  expression  is 


In  general  there  is  no  single,  unique  relationship  between  xc  and  W  . 

3  o 

From  a  fracture  mechanic’s  viewpoint  initial  flaws  of  length  do  not 

alter  the  above  equivalence  expression.  A  typical  example  is  worked  out 
in  Appendix  E  . 
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MULTI-DEGREE-OF- FREEDOM  SYSTEMS 


All  of  the  previously  developed  equati have  expressed  damage  state 
and  rate  parameters  in  terms  of  a  stress  and  stress  cycles  per  time.  These 
same  equations  can  also  be  used  for  multi-degrea-of-freedom  (MDF)  systems 
by  using  the  proper  damage  state  and  rate  parameters.  Reference  (L3)  shows 
that  the  proper  damage  state  and  rate  parameters  are  obtained  by  adding  the 
various  resonant  mode  stresses  and  resonant  frequencies  in  the  mean-square 
sense. 

Consider  the  example  of  a  two-degree-of-freedom  (2DF)  system  whose  stress 
response  is  shown  in  figure  8  . 

fj_  -  center  frequency  of  first  resonant  mode  stress  response  (Ha) 

f2  *  center  frequency  of  second  resonant  mode  stress  response  (Hz) 

2 

S(f)"  stress  power  spectral  denisty  (KSI*/Hz) 

2  2 

o,  *  mean-square  stress  response  of  the  first  resonant  mode  (ESI  ) 

2  2 

<?2  •  mean-square  stress  response  of  the  second  resonant  mode  (KSI  ) 

a„  *  effective  damage  state  stress  (KSI) 

f  .  ■  effective  damage  rate  (Hz) 
eff 

*  J  fb  si  <f>  d  f  (37) 

a 

■J'|d  S2  (f>  d  *  <38> 
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will  take  on  a  value  nearest  the  resonant  mode  having  the  larger 


eff 

stress  power. 


N  ,,  -  f  ,,  x  T 

eff  eff 


(41) 


where  T  ■  test  time 

-  number  of  effective  stress  cycles 


The  above  2DF  case  can  be  extended  to  the  MDF  case  as  follows: 


(42) 


j  -  1,  2,  3  ...  k 

where  j  •  resonant  mode  index 

k  *  total  number  of  resonant  modes 
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RANDOM  VARIABLE  CONSIDERATIONS 

The  results  of  this  study  can  be  characterized  as  being  deterministic. 

That  is,  In  all  cases  the  applied  stresses  and  fatigue  curves  were  treated 
as  deterministic,  not  as  random  variables.  It  is  beyond  the  scope  of  this 
study  to  go  into  the  details  of  how  to  treat  the  random  variable  case. 

It  will  be  pointed  out,  however,  that  methods  have  been  developed 
(e.g.  references  [1],  [2],  [7],  {13  )  for  the  random  variable  case.  The 
values  of  input  acceleration  levels  and  test  time  must  be  converted  to 
stresses  and  stress  cycles  respectively.  The  stress  and  cycles  should  then 
be  treated  as  median  or  average  values.  Standard  deviation  stress  and 
fatigue  curve  values  must  then  be  assigned.  Then  the  random  variable 
expressions  (e.g.  probability  of  failure  versus  cycles)  can  be  used. 
EXAMPLE: 

Given:  63  -  37  Tin-Lead  Solder 
n  »  1 

C5  -  0.292  KSI/G  (SHEAR  STRESS) 

*  7  g  rms  at  accelerated  test  level 

WQ  -  0.025  g2 /Hz 

Quantity  of  solder  joints  being  stressed  •  100 

fa  -  200  Hz 

T2  “  25  minutes  total 

Find:  Average  number  of  cumulative  solder  joint  failures  versus 

test  time. 
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Solution:  r 


2.04  KSI  RMS  shear 


-  0.292  x2  - 

N  ■  12  x  10^  T  cycles 

C  -  6.62  KSI  (See  page  E-12) 

e  -  8.97 


N 

m 


3.6 


x  10 


cycles 


T„  ■  3  minutes 

m 


Choose  the  fatigue  curve  standard  deviation  to  be  102  of  the  median 


value  (i.e.  A/A 
F(T) 


F(T) 


q(T)  -  100  x  F(T) 


(44) 


(45) 


where 

A  •  standard  deviation  of  fatigue  curve  average  value  A 
T  «  test  time  (minutes) 

Tm  -  test  time  for  502  (median)  solder  joint  failures 
q(T)“  average  number  of  cumulative  solder  joint  failures 


as  a  function  of  time. 


Figure  9  is  a  plot  of  q(T)  versus  T  .  About  92%  of  the  solder  joint 
failures  will  occur  after  10  minutes  of  tests. 

If  this  analysis  were  dete minis  tic  (i.e.  A  «  0) ,  all  failures  would 

have  occured  at  m  Tm  ■  3  minutes.  The  scatterband  of  fatigue  curve 

failure  points  results  in  failures  occuring  both  before  and  after  Tffl  . 

Given:  Same  example  as  the  previous  one  involving  63-37  Tin-Lead 
Solder. 

Find:  Equivalent  resonance  dwell  sinusoidal  input  acceleration 

level  x  that  will  produce  the  same  quantity  of  failures 
versus  time  as  shown  in  figure  9. 

Solution:  From  data  on  page  E-12 

-4-  “  2.19 

C 


Using  equation  (98) 

*x'g  -  6.12/qf*  g*s 

JL 

r# 

*s 

-tol 

10 

1.94 

20 

1.37 

30 

1.12 

20 

0.97 

50 


0.87 


FIGURE  9  FAILURE  DISTRIBUTION  FOR  SOLPER 

JOINT  EXAMPLE 
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MULTI  ACCELERATION  FACTOR  CONSIDERATIONS 


In  a  relatively  complex  structural  assembly  such  as  a  typical  electronic 


"black  box"  more  than  one  acceleration  factor  can  be  computed;  one  for 


each  structural  element  that  could  experience  a  fatigue  failure.  Even 


when  only  the  more  highly  stressed  elements  are  given  consideration,  the 


problem  still  exists  of  choosing  the  single  acceleration  factor  to  be 


used  for  testing  the  "black  box." 


At  the  service  environment  the  various  structural  elements  that  could 


fall  will  each  possess  Its  Individual  cumulative  damage  value  (l.e.  potential 


for  failure) .  Due  to  these  cumulative  damage  values  in  conjunction  with 


scatterband  effects  one  structural  element  will  fall  before  the  others. 


If  all  structural  elements  had  the  same  acceleration  factor,  the 


cumulative  damage  and  scatterband  relationship  between  elements  would 


remain  unaltered  at  the  accelerated  test  environment.  This  would  be  the 


case  whether  the  acceleration  factor  Included  non-linear  parameters  in 


Its  estimation  or  not.  The  same  element  would  fall  before  the  others  at 


both  the  accelerated  and  service  environments.  Thus,  scatterband  parameters 


are  not  important  in  estimating  the  acceleration  factor.  Only  the  median 


parameter  values  Influence  the  acceleration  factor  value. 


In  general  all  the  structural  elements  will  not  have  the  same  acceleration 


factor.  Differences  in  damping  linearity,  fatigue  and  crack  growth  rate 


curve  slopes,  and  Initial  flaw  sizes  (n,  6,  6,  &, )  as  examples  will  alter 


the  cumulative  damage  relationship  between  the  structural  elements.  Which 


of  the  various  acceleration  factors  Is  chosen  will  result  In  a  proper 


accelerated  test  for  only  one  class  of  structural  elements.  The  other 
elements  will  either  be  under  or  over-tested.  The  most  conservative 
approach  would  be  to  select  the  largest  acceleration  factor  for  testing 
the  black  box.  The  acceleration  factor  used  could  be  the  average  of 
all  acceleration  factors.  Thus  the  selection  of  a  single  acceleration 
factor  is  considered  to  be  subjective. 
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DECELERATION  FACTOR  CONSIDERATIONS 


It  is  sometimes  of  interest  to  decelerate  a  test  instead  of  accelerating 
one.  An  example  would  involve  a  "black  box"  that  had  been  vibration 
qualified  at  a  relatively  large  input  level  for  a  short  duration  T^. 
It  might  be  desired  to  compute  an  acceptance  test  which  would  produce 
the  same  damage  (i.e.  potential  for  failure)  but  at  a  lower  level  *x2 
for  a  longer  duration  T2 .  In  this  case  *x^  >  xj  and  T2  >  T1 ■  The 
desired  computation  can  proceed  using  previously  developed  equations. 

For  the  fatigue  case  (i.e.  the  initial  flaw  size  a^^  ■  0)  equations  (17) 
or  (111)  can  be  used  as  follows  depending  upon  whether  *x2  or  T2  is  the 
unknown  quantity: 


For  the  fracture  mechanics  case  (i.e.  a ^  >  0)  equations  (105)  and  (106) 

and  program  PL-3  are  applicable.  Equations  (105)  and  (106)  should  be  used 

#• 

if  T2  is  the  unknown  quantity.  PL-3  should  be  used  if  x2  is  the  unknown 

•• 

quantity.  Care  should  be  taken  that  x^  >  x2  and  T2  >  T^.  T^  and  T2  can 
be  in  any  time  units  (e.g.  seconds,  minutes)  so  long  as  T^  has  the  same 
units  as  T2<  As  mentioned  in  a  previous  section  if  the  value  of  a^  is 
chosen  too  large^X  will  be  negative  and  PL-3  won't  execute.  A  negative 
value  of  X  means  that  a^  >  ac,  or  ac^.  That  is,  the  part  will  fracture 
immediately. 


PL-3  deceleration  factor  example: 


Given :  7  07  5-T6  alloy 

a^  •  0.01  inches 

*x^  ■  15  G  rms  (Qual  Test  Level) 

*  1  minute 

T2  *  100  minutes 

Find:  Acceptance  Test  Level  ^  ^or  h  “  0.833,  1 


Using  PL-3  •• 

x2 

n  (G  rms) 

0.833  A. 02 

1  5.61 


As  in  previous  examples  the  deceleration  factor  is  sensitive  to  the 
damping  linearity. 

For  a  complex  "black  box"  where  multi-deceleration  factors  can  exist  the 
choice  of  a  single  factor  is  subjective.  The  most  conservative  approach 
would  be  to  select  the  smallest  deceleration  factor. 
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CONCLUDING  REMARKS 

1.  Expressions  have  been  developed  chat  relate  service  environmental 
level  and  duration  to  accelerated  test  level  and  compressed  duration 
for  the  same  potential  to  do  fatigue  damage.  The  power  law  expression 
commonly  used  in  the  industry  is  shown  to  apply  on  for  those  fatigue 
cases  where  initial  flaws  (i.e.  cracks)  of  length  a^  In  the 
structural  elements  being  stressed  do  not  exist.  Fracture  Mechanics 
effects  (i.e.  where  cracks,  either  actual  or  hypothesized,  already 
exist)  complicate  the  expression;  it  becomes  a  transcendental  function 
whose  solution  is  most  easily  handled  by  the  included  Basic  Language 
computer  program.  Levels  are  in  terms  of  either  stress  or  "black  box" 
vibe  input  acceleration.  Durations  are  in  terms  of  either  number  of 
applied  stress  cycles  or  time. 

2.  The  developed  expressions  apply  to  sine  or  random  vibration  and  thermal 
cycling  for  both  the  low  and  high  cycle  (i.e.  inelastic  and  elastic) 
fatigue  regions.  Linear  and  non-linear  dependence  of  stress  upon 
input  vibration  level  is  included.  Single  and  multi-degree-of-freedom 
systems  are  also  included. 

3.  The  developed  expressions  are  summarized  in  Appendix  G. 

4.  All  equations  are  in  practical  engineering  terms  and  are  expected  to 
be  accurate.  Application  is  straightforward. 

5.  Random  variables  (e.g.  scatterband  fatigue  curve  and  applied  stress)  can 
be  added  to  the  results  of  this  basically  deterministic  study. 
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6.  The  damage  state  and  rate  (i.e.  conditions  of  similitude)  must  be  the 
same  for  both  service  and  accelerated  test  environments. 

7.  Fatigue  damage  is  directly  related  to  the  stress  level  and  number  of 

applied  stress  cycles  in  a  structural  element.  It  is  only  indirectly 

related  to  the  "black  box"  vibe  input  acceleration  level  and  test 

duration.  The  random  response  stress  is  directly  related  to  the  value 

of  the  vibe  input  acceleration  power  spectral  density  W0  in  the 
vicinity  of  the  resonant  frequency.  It  is  only  indirectly  related  to 
rhe  overall  "black  box"  vibe  input  rms  acceleration  level  *x  . 

8.  Two  examples  of  quantitatively  compensating  for  similitude  condition 
violations  are  given. 

9.  The  power  law  relation  applies  only  for  a^  ■  0  .  The  value  of  the 
power  law  exponent  is  l/n6  .  8=2  in  the  low  cycle  fatigue  region 
for  most  structural  materials.  In  the  high  cycle  fatigue  region  8*9 
for  ductile  materials  and  *  20  for  brittle  materials,  n  represents 
the  damping  linearity,  n  «  1  (linear).  0.714  <_  n  <_  2  for  the  cases 
studied. 


A1  <  A  ;  6  <  8 
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0.007  inches 


For  one  example  using  7075-T6  with  - 

I  -  180  KSI  ;  8  -  9.65 

A±  -  111  KSI  ;  6  -  4 


11.  An  initial  flaw  does  not  alter  the  relationship  between  the  sine  and 

random  stress  levels  that  will  propagate  a  crack  of  the  same  size  in  the 
same  time;  namely, 


AS 


(  a1  »  0  )  (49) 


*i 


C 


i 


(50) 

i 


12.  The  random  vibe  transcendental  function  is 

=•-(%)  (+) 

»•  •  • 

where  X  is  a  function  of  both  X2  and 


(51) 


This  equation  cannot  be  normalized  in  closed  form  fashion  as  in  terms 
of  (xj/x^)  and  (T^/^  due  to  the  inherent  non-linearities  in  the 
Fracture  Mechanics  terms. 

13.  A  J-Integral  fracture  mechanics  equation  is  proposed  for  use  in  the 

low  cycle  fatigue  region.  It  has  an  identical  crack  growth  rate  equation 
form  to  the  one  used  in  the  high  cycle  fatigue  region.  Thus,  its 
parameter  values  can  be  substituted  into  the  previously  developed 
high  cycle  equations. 


14.  In  general  there  is  no  single,  unique  relation  between  the  sine  input 

§  * 

vibe  acceleration  level  x^  and  the  random  vibe  input  power  spectral 
density  WQ  .  The  relationship  can  be  altered  as  a  function  material 
ductility,  damping  linearity,  and  fQ/Q  ratio. 

15.  The  developed  equations  also  apply  to  multi-degree-of-freedom  systems 
by  using  the  proper  damage  state  and  rate  parameters.  These  parameters 
are  obtained  by  adding  the  various  mode  resonant  response  stresses  and 
resonant  frequencies  in  the  mean-square  sense.  For  a  2DF  system: 


°T 


1 


fl  <  feff  <  f 2 

feff  *  fl  lf  al  °2 

feff  “  f2  if  °2  »  °1 


16.  The  selection  of  a  single  acceleration  or  deceleration  factor  for  a 
complex  "black  box"  is  considered  to  be  subjective  in  general.  The 
most  conservative  approach  would  be  to  select  the  largest  acceleration 


factor  and  the  smallest  deceleration  factor. 
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SYMBOLS 


Ci  -c5 


crack  half-length 
critical  value  cf  a 


initial  value  of  a 


cross-sectional  area 


material  sine  fatigue  curve  constant 

material  sine  fatigue  curve  constant  with  a^.  >  0 

random  stress  bandwidth 

viscous  damping  coefficient 

constant  of  Paris  crack  growth  rate  curve 

material  random  fatigue  curve  constant 

material  random  fatigue  curve  constant  with  a±  >  0 

constant  of  Dowling-Begley  crack  growth  rate  curve 


constants 


cumulated  damage  at  the  jcn  stress  level 
two-degree-of-freedom 

crack  growth  rate 
modules  of  elasticity 


Error  Function 


frequency 

effective  frequency 
resonant  frequency 


Coulomb  friction  force  value 


cumulative  probability  of  failing  in  time  T 
gravity  acceleration  units 


'  UM 


V  -rf** 


I 


vn 


thickness 


normalized  stress  limit  level 
area  moment  of  inertia 


index 


constant  of  specific  damping  energy  curve 
energy  line  integral 
range  of  J 


critical  value  of  AJ 


total  number  of  resonant  modes,  spring  stiffness 
thousands  of  pounds  per  square  inch 


volume  stress  factor 


stress  intensity  factor  range 

fracture  toughness 

length 


stress  limit  level 


exponent  of  specific  damping  energy  curve 

number  of  stress  cycles 

number  of  stress  cycles  to  failure 

number  of  Independent  opportunities  for  stress  peak  occurrence 
probability 


power  spectral  density 
resonant  amplification  factor 
average  number  of  cumulative  failures 
sinusoidal  stress  amplitude 
sinusoidal  stress  range 


C 


I 


S(f) 


power  spectral  density  as  a  function  of  frequency 
average  time  between  Independent  events 


average  time  between  stress  peaks 


width 


weight 

acceleration  power  spectral  density 


acceleration  level 


dummy  variable 
geometrical  parameter 
response  acceleration 
relative  displacement 
root-mean-square 

fracture  mechanics  dependent  term 


dummy  variable 

fatigue  curve  slope  parameter 


standard  deviation  of  A 


random  rms  stress  value 


constant  of  Paris  crack  growth  rate  curve 
Gamma  Function  with  argument  a 
Incomplete  Gamma  Function  with  argument  a 
applied  strain  amplitude 
fatigue  ductility  coefficient 
material  ultimate  percent  elongation 
applied  strain  range 

average  net  section  plastic  strain  range 
average  net  section  elastic  strain  range 


i 


4 


'If' 


REFERENCES 


1.  R.G.  Lambert,  "Analysis  of  Fatigue  Under  Random  Vibration",  The  Shock 
and  Vibration  Bulletin  46,  Naval  Research  Laboratory,  Washington,  D.C., 
August  1976. 

2.  R.G.  Lambert,  "Mechanical  Reliability  For  Low  Cycle  Fatigue",  Annual 
Reliability  and  Maintainability  Symposium  Proceedings,  January  1978. 

3.  M.  Kitagawa,  "Intergranular  Fatigue  Fracture  of  Chemical  Lead  at  Room 
Temperature",  Cyclic  Stress-Strain  Behavior  Analysis,  Experimentation, 
and  Failure  Prediction,  STP  519,  ASTM. 

4.  Technical  Report  on  Fatigue  Properties  -  SAE  J1099,  February  1975. 

5.  MIL-HDBK-5C ,  Metallic  Materials  and  Elements  for  Aerospace  Vehicle 
Structures,  September  1976,  Superintendent  of  Documents ,  Government 
Printing  Office,  Washington,  D.C. 

6.  Part  1,  Damage  Tolerant  Design  Handbook,  MCXC-HB-01,  Metal  and 
Ceramics  Information  Center,  Battelle  Columbus  Laboratory,  Columbus, 
Ohio. 

7.  R.G.  Lambert,  "Fracture  Mechanics  Applied  to  Step-Stress  Fatigue  Under 
Sine/Random  Vibration",  The  Shock  and  Vibration  Bulletin  48,  Part  3, 
Naval  Research  Laboratory,  Washington,  D.C.,  September  1978. 

8.  N.E.  Dowling  and  J.A.  Begley,  "Fatigue  Crack  Growth  During  Gross 
Plasticity  and  the  J-Integral,"  ASTM  STP  590,  1976,  p.  83. 

9.  D.  Broek,  "Elementary  Engineering  Fracture  Mechanics,”  Sijthoff  & 
Hoordhoff,  Alphen  aan  den  Rijn  -  The  Netherlands,  1978. 

10.  D.F.  Mowbray,  "Derivation  of  a  Low-Cycle  Fatigue  Relationship 
Employing  the  J-Integral  and  Associated  Crack  Growth  Law,"  ASTM  STP 
601,  1976,  p.33. 

11.  D.F.  Mowbray,  "Study  of  Fatigue  Crack  Growth  at  High  Growth  Rates," 
Materials  and  Processes  Laboratory,  General  Electric  Company, 
Schenectady,  New  York,  July  1977  (unpublished  work). 

12.  S.  Rolfe  and  J.  Barsom,  "Fracture  and  Fatigue  Control  In  Structures— 
Applications  of  Fracture  Mechanics",  Prentice-Hall,  Inc.,  1977. 

13.  R.G.  Lambert,  "Fatigue  Analysis  of  Multl-Degree-of-Freedom  Systems 
Under  Random  Vibration",  The  Shock  and  Vibration  Bulletin  47,  Naval 
Research  Laboratory,  Washington,  D.C.,  September  1977. 


57 


14.  M.  Gertel,  "Specification  of  Laboratory  Tests"  pp  24-9  to  24-10, 
Volume  2,  Shock  and  Vibration  Handbook,  edited  by  Harris  and  Crede, 
McGraw-Hill  Book  Company,  1961. 

15.  C.  Crede  and  J.  Ruzicka,  "Theory  of  Vibration  Isolation",  pp  30-7  to 
30-8,  Volume  2,  Shock  and  Vibration  Handbook,  edited  by  Harris  and 
Crede,  McGraw-Hill  Book  Company,  1961. 

16.  B.  Lazan  and  L.  Goodman,  "Material  and  Interface  Damping" 

pp  36-33  to  36-35,  Volume  2,  Shock  and  Vibration  Handbook,  edited 
by  Harris  and  Crede,  McGraw-Hill  Book  Company,  1961. 


APPENDIX  A 


E» 
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SIMILITUDE  COMPENSATION  EXAMPLES 

The  equations  developed  in  this  study  to  calculate  test  levels  to  cumulate 
the  same  fatigue  damage  per  the  desired  time  compression  ratio  are  accurate 
only  if  the  conditions  of  similitude  (i.e.  damage  state  and  rate)  are  the 
same  in  both  environments.  There  are  many  common  ways  for  such  conditions 
to  be  violated.  Only  one  such  violated  condition  and  a  technique  to 
compensate  for  the  violation  will  be  treated  in  this  section. 

Reference  [1]  shows  that  most  fatigue  damage  during  random  vibration  is 
caused  by  stress  peaks  between  2a  to  5o  (where  o  ■  rms  stress  level) 
for  ductile  materials  and  between  3 a  to  6cr  for  brittle  materials.  If 
the  peak  stresses  are  limited  from  exceeding  some  limit  level  L  (KSI) ,  the 
fatigue  life  will  be  extended  from  that  for  no  stress  limiting  (i.e.  L  ■  ») . 

Electronic  equipment  mounted  in  rocket  propelled  spacecraft  or  jet  aircraft 
experience  "unlimited"  stress  peaks.  Note:  for  a  frequency  band-limited 
process  the  very  high  (i.e.  >  6c)  peaks  do  not  occur  very  often  even  when 

limiting  is  not  present.  For  a  typical  stress  response  process  that  is 
band-limited  to  500  Hz:  6c  peaks  occur  about  every  36  hours  and  7c  peaks 
occur  about  every  3  years  on  the  average.  Electrodynamic  shaker  systems  also 
produce  high  peaks  in  a  gaussian  fashion  due  to  peak  restoration  by  the  shaker 
transfer  function  as  long  as  the  noise  generator  voltage  is  hard  limited 
(i.e.  clipped)  no  lower  than  3c.  Such  peaks  may  not  occur  in  a  test  of 
relatively  short  duration.  See  Appendix  H. 

Stress  limiting  will  occur  if  the  noise  generator  voltage  of  an  electro- 
dynamic  shaker  system  is  clipped  below  3c  ,  if  the  shaker  system  (mechanical 
or  otherwise)  is  incapable  of  producing  5  or  6c  peaks  conceivably  due  to 
Brinelling  of  metal  surfaces  or  compressibility  of  fluids,  or  if  the 
structural  element's  motion  being  vibrated  is  snubbed  by  design. 

A-l 


The  effective  limit  level  is  L/a  . 


Define  H 


2+6 


(52) 


where 


L  - 


a 


e  - 


Nf  - 


N*  “ 


limit  stress  level  (KSI) 

rms  stress  level  (KSI  BMS) 

material  fatigue  curve  slope  parameter 

fatigue  life  (i.e.  cycles  to  failure)  for  no 

limiting  (i.e.  L  ■  ®) 

extended  fatigue  life  with  stress  limiting 


From  reference  [1] 


r  (q) 


T  (a  ,  H)  + 


r  (a) 

Y  (a) 


Gamma  Function 
Incomplete  Gamma  Function 


(53) 


(54) 


The  following  table  shows  the  factor  of  fatigue  life  extension  for  various 

limit  level  on  Copper  wire. 

TABLE  VII  LIMIT  LEVEL  FATIGUE  LIFE  EXTENSION 


Limit 

Level 

(L/a) 

»£l  /  »f 

00 

1 

4 

1.08 

3 

1.86 

2 

13.0 

A-2 


For  the  purposes  of  this  analysis  assume  that  no  limiting  occurs  at  service 
vibration  levels  but  that  L/c  *  3  at  the  accelerated  test  level.  The 
compensation  technique  is  to  extend  the  desired  test  time  by  the  factor 
Nf^/Nf  .  In  the  above  example  T2  would  be  multiplied  by  1.86. 

A  second  type  of  similitude  condition  violation  compensation  method  will 
be  treated. 

Given:  The  similitude  violation  is  due  to  a  difference  in  stress  spectra 

between  service  and  test  environments.  The  stress  vibration  system 
is  2DF.  The  calculated  acceleration  factor  (j^/x^)  is  3  for  the 
desired  time  compression  factor  O^/Tj)  using  equation  (111).  The 
spectra  parameter  values  are  given  below  and  in  Figures  10  and  11: 


PARAMETER 

ENVIRONMENT  | 

SERVICE 

TEST 

oa  (ksi) 

8 

24 

ob  (ksi) 

16 

20 

fa  (Hz) 

150 

150 

*b  (Hz) 

375 

375 

Op  (ksi) 

17.9 

31.2 

fef f  (Hz) 

342 

264 

Find:  The  appropriate  test  compensation  factors  such  that  oT 


TEST 


3  0, 


L SERVICE 


and  f 


eff, 


-  f 


TEST 


eff 


SERVICE 


Solution:  It  can  be  seen  that  the  resonant  frequencies  are  the  same  at 

both  environments.  However,  did  not  increase  from  16  to  48  ksi 

as  desired.  This  would  cause  an  inappropriate  test  damage  state 

and  rate.  oT  should  be  3  x  17.9  ■  53.7  ksi.  Therefore,  x« 
iTEST  2 


A-3 


needs  to  be  increased  by  an  additional  factor  of  53.7/31.2  ■ 
1.72.  The  test  duration  needs  to  be  increased  from  its  computed 
compressed  value  by  a  factor  of  342/264  ■  1.3. 

Alternative  Solution: 

In  some  cases  it  is  possible  to  increase  o^  at  the  test 
environment  from  20  ksi  to  the  desired  value  of  48  ksi  by 
the  use  of  shaker  equalization  techniques.  The  spectra  obtained 
is  the  desired  test  spectra.  No  further  compensation  to  overall 
rms  level  or  test  duration  is  required.  This  is  the  most  direct 
and  preferred  solution. 


FIGURE  10 


SERVICE  LEVEL  STRESS  SPECTRUM 


FiG-URE  II  TEST  LEVEL  STRESS  SPECTRUM 
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APPENDIX  B 


DAMPING  TEEM  DERIVATION 


The  power  lew  exponent  is  called  the  damping  term  n  when  related  the 
stress  at  an  internal  structural  element  due  to  the  "black  box"  input 
vibe  acceleration  level: 

SINE:  ~  -  S  -  c2  S  (55) 

m  nR 

RANDOM:  a  -  C4  x^  (56) 


The  relationship  between  stress  and  acceleration  is  linear  for  n  »  1  . 

Non-linearities  in  this  relationship  can  arise  if  any  of  the  elements  of 
the  idealized  mass-sprlng-dashpot  system  becomes  amplitude  sensitive 
(i.e.  non-linear).  Several  such  cases  will  be  evaluated.  The  technique 
used  in  reference  P-4j  for  developing  the  n  determination  will  be  used  here. 

CASE  A  (Q  -  constant:  sinusoidal) 

y  z  • 

S  - 
A  ■ 

x  •• 

At  resonance: 

*•  •• 

my  *  SA  ;  mQx  ”  SA 

s  -  A  *x 

A 

ns"  1  0 


y  -  x 
stress 

spring  cross-sectional  area 
(57) 


B-l 


E  (Coulomb  Friction:  sinusoidal) 


For  Q  >  3  ;  Reference  I 15] 


Q  * 


2k* 

£  4ir2  f, 


(58) 


my  ■  mQx  »  SA 


SA  SAF.  4ir2  f  2 

“Q  2  m  k  x  ” 

’  2  0 


2km 

T 


A  Fj  Air 


~2 

x 


(59) 


C 


(internal  sinusoidal  stress-strain  hysteresis) 
From  reference  {14J,  [.163 


z  -  y  -  x 
S  -  C2  2 


Kv  ir  S2 
E  JeSn 


n  *  2.4  for  most  structural  materials  (6i 

stressed  below  0.8  of  fatigue  strength 


mQx 


SA  \  n 


8  for  higher  stresses 


CASE  D  (internal  random  stress-strain  hysteresis) 


Same  figure  as  for  CASE  C 


••  •# 

z  =  y 
rms  J  rms 


m  ir^  fn  Kv 

A  2EJo  (fb-fa) 


1/n 


2/n 


-  -  -  ~  -  0.833 
n  2.4 


'R 


0.833 


for  structural  materials 


For  viscoelastic  adhesives  stressed  in  shear  £14} 

"»  ‘  dr  '  °'78‘  <? 

Figure  12  shows  the  non-linear  relationship  between  a  and  z  for 
C4  »  1.0  KSI/G  . 


I 

i 


I 


1 


CASE  E  (Coulomb  Friction:  random) 
Same  figure  as  for  CASE  B 


“rms 


yrms 


m^rms  * 
Frms 


'rms 


Q  * 


*f%~{ ~o  Q  ' 
A  / ff-  fp  £ 

7  2  <*b  ‘  V 


2  k  x 


ff  ,™2  £„2 


1  / 

_ 


‘rms 


'•rms 


aJ  mo  <*b  * 


fa> 


g  A 
m 


Ff  <fb  -  fa>  fn 


1.5 


nR  "  1,5 


CASE  F  (Viscoelastic  Materials:  sinusoidal) 

Same  figure  as  for  CASE  C 

For  viscoeleastic  materials  stressed  in  shear  [14]  n  ■  2.55. 
equation  (62) 


Tic  *  - 1 -  »  . 

s  n  -  1  2.55  -  1 

ns  -  0.645  Q 


(65) 


From 


B-4 


('ksl  rms) 


B-5 


APPENDIX  C 


FRACTURE  MECHANICS  FATIGUE  CURVE 
PARAMETER  APPROXIMATIONS 


(6  -  2)  cfl  Y°  AS' 


Ct) 


e  -  2 
2 


from  reference  [7]  (66) 


(0  -  2) 


S  »f  2®  (  >i) 


6-2 

2 


Define 


A^  *  sinusoidal  curve  "y  -  intercept"  parameter  with  an  initial 
flaw  of  length  ai 


*  1 


(6  -  2)  cQ  Y"  2' 


(-i) 


6  -  2 
2 


»-(i) 
•  •(+) 


e  „ 


(6  -  2)  cQ  Ye  Nf2 


fc) 


e  -  2 
2 


-fa. 


f «.  fc.”  v  ■ 


Mm 


Define 


C 


i 


random  fatigue  curve  "y  -  intercept"  parameter  with  an  initial 
flaw  of  length  a^^ 


C 


i 


ksi 


(70) 


(71) 


7075-T6 


1/9  -  0.25  ;  C  ■  6  x  10~9  in/cycle  ;  Y  ■  1 

AKC  *  20  KSI  IN.' 


AS  a  F  1.06125 

2  L  &i  Nj 


T  l/e 


x  ioc 


ksi 


"y  -  intercept" 


I  1.06125  x  106 

A*  “  L — * — J 


AS/ 2  « 

0.25 


Ai  J  Nf 


ksi 


(72) 


TABLE  VIII  A±/\  RATIO  VALUES 


(inches) 

A^  (ksi) 

C±  (ksi) 

V‘i 

.007 

U1 

49.3 

2.25 

.050 

67.9 

30.2 

2.25 

.100 

57.1 

25.4 

2.25 

A  •  180  ksi 
"C  -  80  ksi 


-  -  2.25 


C-2 
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APPENDIX  D 


FRACTURE  MECHANICS  ACCELERATED  SINUSOIDAL  TEST  STRESS  LEVEL  DERIVATION 


From  references  [1]  and  [7] 


E1  “  (e  -  2)  c0  ASj"  Y8 


6  -  2  6-2 


ft)  ‘  -  ft) 


(73) 


Nf2  "  (9  -  2)  cQ  AS2b  Y° 


0-2  6  -  2 


ft)  ‘  -ft) 


(74) 


%  ■  [alj]  ;  S  ■  [St]  <75> 

where 

■  number  of  cycles  Co  failure  ac  service  level  1 
Nf^  ■  number  of  cycles  to  failure  at  accelerated  test  level  2 
6  »  constant  of  crack  growth  rate  curve 

Y  “  geometrical  parameter 
cQ  ■  constant  of  crack  growth  rate  curve 


D-l 


a  “  initial  crack  length 

ac  -  critical  crack  length 

AKc  *  fracture  toughness 

AS  «  applied  stress  range 

Define  -  number  of  applied  stress  cycles  at  level  1 

N0  *  number  of  applied  stress  cycles  at  level  2 

n2 


For  equal  damage 


(78) 


!i  A i\ 

Ni  “  vas2  / 


(■) 


where 

X 


correction  factor;  X  corrects  for  the  dependence  upon  ai  ,  Y 


AS ^  >  AS £  9  AKC  9  9 


9-2 

&Y 


(*) 


9  -  2 
2 


(79) 


Substituting  from  previous  equations 

*i  Y2  ASX2 


AK„ 


) 


9  ~  2 
2 


at  Y2  AS22 


9-2 

2 


AK, 


(80) 


X  >  1  for  AS2  >  ASX 


Let  h 


(*?) 


8-2 


8-2 


L  -  1  -  h  AS, 


(81) 


(82) 
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The  above  two  equations  are  transcendental  functions  and  must  be  solved 
accordingly. 

AS 2  >  AS1 

Nx  >  N2 

X  >1 

A  negative  value  of  X  indicates  that  a^  >  ac^  or  a^  >  a^  .  A  zero 
value  of  X  indicates  that  ■  ac^  .  Such  situations  are  unrealistic 


for  this  analysis.  Fracture  would  occur  during  the  application  of  the  first 
stress  cycle  of  either  AS.  or  AS,  . 


X  will  be  negative  if  the  selected  value  for  is  larger  than 


AS- 


max 


'  ^c2  ] 

_  "1  ?2  J 


1/2 


(85) 


Thus  the  value  of  AS2  that  will  satisfy  the  transcendental  equation  is 


<_  as2 


max 
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APPENDIX  E 


SINE-RANDOM  EQUIVALENCE  DERIVATION 


It  is  sometimes  desired  to  determine  the  sinusoidal  input  acceleration 

*x_  (g's)  at  a  structural  resonant  frequency  f  (Hz)  that  will  cumulate 
s  n 

the  same  fatigue  damage  as  a  wide  band  random  vibe  input  of  power  spectral 
density  W  (gz/Hz).  x  is  sometimes  [17)  referred  to  as  the  equivalent  sine 


i\ 


input.  It  'will  also  be  determined  if  a  single,  unigue  x  -W_  relationship 
exists  for  all  structural  elements  in  all  "black  boxes".  The  assumptions 
made  are  that  the  structural  elements  being  stressed  at  resonance  can  be 
characterized  as  single-degree-of-freedom  systems  and  that  the  duration 
of  both  sine  and  random  vibe  tests  are  the  same.  Patigue  and  Fracture 
Mechanics  effects  are  considered. 


Consider  the  following  idealized  single-degree-of-freedom  system: 


71  _J 


Sj<f) 


(g2/Hz) 


z  *  strech  in  spring  ■  y  -  x 


fn  -  _L  A  fP 

n  2ir  y  m 


fa  <  f_  <  f. 
a  n  b 


f  (Hz) 


TrzftHi 


r 


DO 

J  Sj.(£)df  - 


VV£a> 


x  *  random  vibe  input  acceleration  (g's  rms) 

2 

WQ  *  Power  Spectral  Density,  PSD  (g  /Hz) 
fa,f^  *  frequency  limits  (Hz) 


The  new  narrow-band  response  is: 


z  *  y 


9.8  y. 


f  fnWoQ 


i-X-  Q 

XX 


(inches  rms) 


(g’s  rms) 


at  resonance 


12.28  /  „  I  ,  . 

— JTJ  a!  Wq  Q  (inches  rms) 


a  -  rms  stress  -  C,  y _ _ 

3  rms 

C3  “  configuration  constant  (RSI  RMS /INCH  RMS) 


0  “  Co  /  1 r 


'WfW  Wc 


(RSI  RMS) 


It  should  be  noted  that  0  is  directly  related  to  the  value  of  Wq  in  the 

vicinity  of  the  resonant  frequency  fn  (i.e.  between  the  half-power  points 

•  • 

of  the  response  curve),  a  is  only  indirectly  related  to  x^^. 

**rms  "  V  W£.>  '  <«» 

•  • 

The  value  of  can  be  changed  by  altering  the  value  of  WQ  outside  the 

vicinity  of  fQ.  The  magnitude  of  0  will  not  change  significantly.  Thus 


V*  , 


there  is  not  a  unique  relationship  between 
a  and  the  value  of  WQ  in  the  vicinity  of  fn. 


and  a.  There  is  between 


For  sine  resonance  dwell: 


S 


AS 


T 


(KSI) 


S  ■  stress  amplitude 

AS  •  peak-peak  stress  range 
•  • 

xg  -  sine  vibe  input  acceleration  (g  s) 

S  -  z 

■  configuration  constant  (KSI/INCH) 


An  example  will  be  worked  out  to  illustrate  the  use  of  the  above  expression. 
The  system  considered  will  be  that  of  a  massless  beam  of  rectangular  cross 
section  with  a  concentrated  mass  load. 


W  »  weight  (lbs) 

2 

m  ■  mass  (lb  -  sec  /in) 

2 

g  ■  acceleration  of  gravity  ■  386  in/sec 
E  ■  modules  of  elasticity  (lb/in^) 

I  -  area  moment  of  inertia  of  beam  cross  section  (inches^) 


For  a  rectangular  cross-section 


V:  >■»*• 


VJ  -  0.2  g^ /Hz 

o  ^ 


nns 


19.9  g  rms 


nas 


19.8  g  rms 


y  “  0.0777  inches  rms 
rms 


a  -  C3  y 


rms 


(389) (0.0777)  -  30.2  KSI  rms 


Also 


12.28  0 


1/2 


rms 


£»1,5«b  -  £.)1/2 


rms 


Define 


a/2 


C4  „  Cl  ^,12.28}qJ 

^5(fb  -  fa)1/2 


KSI  RMS 
g  BMS 


C4  -  1.517 


Cl 


.  l-5Eh 


a  -  1.517  x  (KSI  RMS) 
rms 


For  Che  sine  resonance  dwell: 


-  f.  2 


x„  -  *s_Ja_  g's 
8  9.8 


s 


e. 


»  r.O-r 


CiQ  9.8  x. 


Define 


c2 


:x  Q  9.8 

— n 


38.1  MI 

g 


S  *  C2  xs 

s  -  38.1  ;; 

For  xs  -  lg  ;  S  -  38.1  KSI 


In  summary  for  che  example  given  where  Wq 


0.2  g2 /Hz  ;  xs  -  lg 


a  -  30.2  KSI  RMS 


S  ■  38.1  KSI  vector 


From  reference  12] 

/  80  \  9,65 
^30.2  J 

■■  ■(+)'  ■ 

(ft)  5'65 

1.15  x  104  cycles  (95) 


3.22  x  106  cycles  (96) 


where  and  Ng  are  cycles  Co  failure  during  the  random  and  sine  tests 


respectively. 


E-7 


rr- 


ss 


Thus  the  value  of 


and  WQ  chosen  for  the  example  are  not  equivalent. 


The  general 


equivalence  case  will  be  continued  as  follows: 


S  -  £1  «  A  No  "1/6 

2  s 

s  -  c2 

NS  “  fn  Ts 

o  ■  c  nr'1/S 

o  “  C*^ 


0 

zrms 


C 

4 


C:  (12.28)  01/2 
f„3/2(t-f  )1/2 

R  D  ft 


Cx  Q  9.8 


T)  -  0.714 

nR  -  0.833 


non-linear  damping  parameter  for 
internal  stress-strain  hysteresis 
damping 


For  the  sine  and  random  tests  to  be  of  the  same  duration,  Tg  ■  TR 


or  N$  -  Nr 


E-8 


r 

fcJ 

% 

'  y 


s  . 

A 


C2  S 

A 

*» 

C4  XR  R 

C 

'  W/F  (^l) 

O  -  ['.  <w]  V2 


Using  previously  defined  expressions  It  can  be  shown  chat  the  desired 
equivalency  expression  Is 


••  ns 

x„  «  (1.25) 


G)i/F  [a-  J 


1  -  nR 


The  above  expression  Is  a  function  of  ng  and  nR.  In  general  ng  and 

nR  values  are  not  restricted  to  any  particular  value.  Thus  from  that 

standpoint  it  can  be  concluded  that  for  structures  with  non-llnear 

— 

damping  there  is  no  single,  unigue  relationship  between  xg  and  WQ. 

Consider  the  linear  damping  case  where  ng  ■  nR  **  1.  The  equivalency 

equation  becomes: 


W„ 


From  [1] 


(1.25) 


•f  *  V7  [r(^)] 


.  .  >  C'X k**" 


I 


*1 

V*  , 


where  6  ■  slope  parameter  of  Che  material's  sinusoidal  fatigue  curve 

-  takes  on  values  of  about  2  for  ductile  materials  and  about  3  for 
brittle  materials. 

TABLE  IX  EQUIVALENCY  RATIO  VALUES 


MATERIALS 

Vw01/2 

7075-T6  (ductile) 

A231-B  (brittle) 

2.81 

3.98 

Thus  it  can  be  concluded  that,  if  a  "black  box"  contains  ductile  and 
brittle  materials,  there  is  no  single,  unique  xs/W0  relationship. 

It  can  also  be  seen  that  the  equivalency  expression  is  a  function  of 
fn/Q.  In  general  f^/Q  is  not  a  constant.  This  fact  also  illustrates  a 
lack  of  a  unique  equivalency  relationship. 


From  Fracture  Mechanics  Considerations  I 7J 

A^  »  sinusoidal  fatigue  curve  "y  -  intercept"  parameter  with  an  initial 
flaw  of  length  a^ 

Ci  -  random  fatigue  curve  "y  -  intercept"  parameter  with  an  initial 
flaw  of  length  a^ 


E-10 


h  _A_  Xi  +  1 

^  C  *  C±  +  C 

Initial  flaws  do  aoc  alter  the  j!g/W0  relationship. 

For  63  -  37  Tin-Lead  Solder  in  the  high  cycle  fatigue  region  at  room 
temperature : 

Reversed  bending: 

A  -  15.3  KSI  ;  c"  ■  6.74  KSI 

6  -  9.85  ;  A/C  -  2.27 

Reversed  shear: 

A  -  14.5  KSI  ;  C  -  6.62  KSI 

8  -  8.97  ;  A/C  -  2.19 

For  Copper  Wire: 

A  -  81.9  KSI  ;  C  «  36.9  KSI 

8  -  9.28  ;  A/C  -  2.22 

The  above  solder  and  copper  can  be  considered  as  ductile.  If  workmanship 
defects  are  treated  as  having  Initial  flaws  of  some  abltrary  lengths,  then 

Ai  <  A  and  <  C  . 


E-ll 


I 


However 


The  equivalency  equation  becomes 


(100) 


In  general  ^2.  £  constant.  The  structural  damping  (hence,  Q)  is 

composed  of  several  damping  mechanisms  (e.g.  friction,  internal  stress- 
strain  hysteresis,  air,  viscoelastic).  A  single  equivalence  still  does 
not  exist. 


Q  *  20  is  a  typical  value  for  an  MLB  without  any  special  added 

viscoelastic  damping  treatment  for  fn’s  ranging  from  100  -  300  Hz. 


For  If  -  0.1  gVHz 

o 


Xg  -  0.197 


g's  (See  Figure  13) 


<3=20 

Wc  =  O-  I  9  /hy 


(101) 


<?*) 


100  (50  ZOO  Z50  300 

^/n  C 

FI CUBE  L3  EQUIVALENT  xg 
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APPENDIX  F 

DERIVATION  OF  ACCELERATED  RANDOM  VIBE  TEST  ACCELERATION  LEVEL 

Refer  to  Appendix  E  for  the  definitions  and  units  of  the  stress  terns 

#»  *  • 

,  02  ,  AS^  ,  and  AS 2  ,  the  acceleration  terms  x^  and  X2  and  the 
stress  per  ms  acceleration  tern  C^  . 

For  a  given  structural  assembly  where  the  shape  of  the  input  power  spectral 
density  S«  (f)  and  the  structural  resonant  frequency  is  the  same  at  both 
levels  1  and  2: 

..  *  m* 

al  “  C4  xi  i  °2  '  C4  x2 


The  equivalent  sine  stress  levels  are: 
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I 


(108) 


X  will  be  negative  if  the  selected  value  for  x2  is  larger  than 


max 


AKC 
ai  c5 


(110) 
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It  can  be  seen  that  the  damping  term  n  alters  the  effective  value  of  the 
crack  growth  rate  parameter  0  .  That  is,  n0  will  be  greater  or  less  than 
0  if  n  is  greater  or  less  than  unity. 


APPENDIX  G 


SUMMARY  OF  DEVELOPED  EQUATIONS 


FATIGUE  (  ai  -  0  ) 


— 

x2 

vr- 

*1 


(SINE  OR  RANDOM)  (111) 


(SINE) 


(112) 


fa  . 

al 


(RANDOM)  (113) 


N  -  fT 


(114) 


AC£ 

A*1 


(CYCLIC)  (115) 


FRACTURE  MECHANICS  (  >  0  ) 
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APPENDIX  H 


TIME  BETWEEN  STRESS  PEAKS 

All  of  che  high  level  random  stress  peaks  may  not  occur  in  a  test  of 
finite  duration.  The  primary  reason  Is  that  such  high  peaks  are  low 
probability  of  occurrence  events  and  don't  occur  very  often.  This 
section  develops  the  approximate  relationship  between  the  average  time 
between  occurrences  T  of  stress  peaks  L.  L  will  be  expressed  in  multiples 
a  of  the  rms  stress  level  a.  (i.e.  a  ■  L/a) . 

The  envelope  of  positive  stress  has  a  Rayleigh  probability  density  function. 
For  values  of  stress  greater  than  about  2. So  the  stress  envelope  approximates 
the  stress  peak  values.  The  percent  of  time  (i.e.  probability)  that  the 
stress  envelope  is  equal  to  or  greater  than  a  can  be  expressed  as 

p  -  e"°2/2  (124) 

o  -  L/o  (125) 


For  a  random  stress  process  of  bandwidth  B  (Hz)  the  average  time  between 
Independent  events  is  approximated  as 

Cj  *  — —  seconds  (126) 

1  B 


The  number  of  Independent  opportunities  for  occurrence  is 


N 


opp 


T 


Ci  N0PP 


H 

p 


T 


seconds 


(127) 


(128) 


-49* 


T 


The  following  conversion  f&ccors  can  be  used  with  T  : 

60  seconds /minute 
3600  seconds /hour 
720  hours /month 
8640  hours /year 


For  B  -  500  Hz  : 


a 

T 

3 

0.18 

sec 

3.5 

0.91 

sec 

4 

6.0 

sec 

4.5 

50.0 

sec 

5 

9.0 

minutes 

5.5 

2.0 

hours 

6 

36 

hours  (1.5  days) 

6.5 

1.2 

months 

7 

2.8 

years 

Equation  (128)  can  be  rearranged  as  follows: 
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APPENDIX  I 


SUMMARY  REPORT 


This  appendix  summarizes  Che  more  often  used  sections  of  the  entire  study. 
The  equation  and  figure  numbers  are  independent  of  those  in  the  main  body 
of  the  report. 


It  2lS  typically  of  great  interest  and  practical  importance  to  accelerate 
a  Gaussian  random  vibration  qualification  or  acceptance  test  of  an 
electronic  "black  box"  in  the  laboratory  from  the  actual  service  conditions. 
The  test  duration  is  compressed  by  a  relatively  large  factor  (e.g.  1000) 
with  an  attendant  increase  in  the  applied  vibration  level. 

Mathematical  relationships  have  been  developed  [1]  which  determine  the 
proper  increase  in  the  vibration  input  root-mean-square  (rms)  level  to  the 
electronic  "black  box"  for  the  desired  test  duration  (i.e.  time)  compression 
factor  such  that  the  accumulated  fatigue  damage  is  the  same  for  both  the  test 
and  service  environments.  Cumulative  fatigue  damage  does  not  necessarily  mean 
structural  fracture  or  failure.  It  means  that  useful  life  is  being  consumed 
and  indicates  the  potential  for  failure.  The  failure  potential  must  be  the 
same  at  both  environments. 


This  paper  shows  the  criteria  for  selecting  the  form  of  the  input  vibration 
level-duration  relationship  and  assigning  values  to  the  parameters.  Derivations 
are  shown  in  reference  [1]. 


r  fv  V"  ’ 


CONDITIONS  OF  SIMILITUDE 


Certain  conditions  of  similitude  must  be  imposed  upon  the  service  and 
laboratory  accelerated  test  environments  if  the  developed  mathematical 
relationships  are  to  be  appropriately  and  accurately  applied.  The 
fundamental  hypothesis  is  that  the  damage  states  and  damage  rates  must 
be  the  same  far  both  environments.  Specifically  the  states  of  stress 
(torsion,  bending,  axial),  the  corresponding  fatigue  strengths,  the 
resonant  mode  shapes,  the  internal  response  stress  spectrum  shapes,  the 
stress  peak  distribution,  and  the  type  and  location  of  failure  mechanisms 
must  be  the  same  for  both  environments. 

Extreme  temperature,  humidity  or  corrosive  element  differences  between 
the  service  and  test  environments  may  result  in  similitude  violations,  if 
such  differences  are  sufficient  to  alter  the  material's  fatigue  strength 
parameters.  Threshold  sensitive  or  other  non-linear  response  effects  in 
general  tend  to  violate  conditions  of  similitude.  In  some  cases  lack  of 
similitude  can  be  quantitatively  compensated  for.  Several  examples  are 
Included . 

The  condition  that  the  shape  of  the  vibration  acceleration  input  spectra 
cr  the  overall  acceleration  rms  levels  must  be  the  same  for  both  environments 
has  purposely  been  emitted  from  the  previously  listed  conditions.  This  is 
because  the  fatigue  damage  state  and  rate  are  only  indirectly  related  to  the 
input  acceleration  spectrum.  They  are  directly  related  to  the  response  stress 
spectrum  at  the  location ■ where  damage  is  accumulating.  The  response  stress 
3pectrum  is  related  to  the  vibration  acceleration  power  spectral  density  value 
in  the  vicinity  of  resonances. 


1-2 


■w  **« 


t 


*1 

.  •* 

t 


CUMULATIVE  FATIGUE 

Black  box  structural  elements  (e.g.  solder  joints,  vires,  device  leads, 
support  structures)  that  are  subjected  to  random  vibration  loads  will 
always  cumulate  fatigue  damage.  Such  fatigue  damage  can  range  In  value 
from  very  little  to  very  large.  It  Is  never  zero.  Structure  element 
fracture  occurs  when  the  cumulated  fatigue  damage  becomes  large  and  is 
defined  mathematically  by  the  material's  fatigue  curve. 


It  can  be  shown  that  curves  of  equal  damage  have  the  same  slope  as  the 
material's  fatigue  curve.  Fatigue  curve  parameter  values  are  readily 
available  from  many  published  sources  (e.g.  MIL-HDBK-5C ,  SAE  J1099) . 

Therefore  it  is  useful  to  work  with  fatigue  curve  parameters  for  determining 
test  acceleration  factors  whether  or  not  large  fatigue  damage  is  accumulated. 

For  sinusoidal  stressing  the  material's  fatigue  curve  is  expressed  as: 

S  >  A  Ns‘1/6  (1) 


where 


Stress  amplitude  (ksi) 

number  of  sinusoidal  stress  cycles  to  failure 
y-intercept  on  log  -  log  plot  for  Ng  ■  1; 
true  ultimate  stress  (ksi) 
slope  parameter 


For  Gaussian  random  stressing  the  material's  fatigue  curve  is  expressed  as 


C  N  ”1/B 
a 


VV*-  •  ti 


where 


a  • 


rms  stress  value  (ksl) 


N_  -  median  cycles  to  failure 
m 

C  «  y-intercept  on  log  -  log  plot  (ksi) 
3  -  slope  parameter  (see  Table  I) 


The  random  fatigue  curve  parameters  can  be  obtained  from  the  sinusoidal 
fatigue  curve.  It  should  be  noted  chat  3*9  for  ductile  materials  and 
S  a  20  for  brittle  materials  regardless  of  the  material's  ultimate  strength. 


Equation  (2)  shows  that  the  fatigue  process  is  directly  related  to  the  rms 
stress  level  o  of  the  structural  element  (e.g.  solder  joint)  inside  the 
electronic  black  box  being  stressed  and  the  median  number  of  applied  stress 
cycles  Nm  .  Both  a  and  are  most  frequently  estimated  by  determining  the 
resonant  response  of  all  the  structural  members  of  the  black  box  and  ad j accent 
structures. 


TABLE  I  TYPICAL  3  VALUES 


MATERIAL 

6 

Copper  Wire 

9.23 

Aluminum  Alloy: 

6061-T6 

7075-T6 

8.92 

9.65 

Soft  Solder 
(63-37  Tin-Lead) 

9.85 

4340  (BHN  243) 

10.5 

4340  (BHN  250) 

13.2 

AZ313  Magnesium  Alloy 

22.4 
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FATIGUE-TEST  RELATIONSHIPS 

Ic  was  previously  shown  chat  the  fatigue  process  was  directly  related  to  the 
rms  stress  level  a  and  the  median  number  of  applied  stress  cycles  Nm.  For 
purposes  of  test  it  is  of  more  interest  to  express  cumulative  damage  in  terms 
of  vibration  input  acceleration  rms  level  x  to  the  electronic  black  box  and 
test  duration  T.  c  can  be  related  to  jc  as  follows: 


C4  *  constant  (ksi/g  rms) 
n  *  damping  parameter 

Table  II  shows  typical  n  values.  Figure  1  shows  typical  a-x  relationships 
for  various  n  values  with  C4  -  1.0  ksi/g.  The  value  of  n  is  best  determined 
empirically  because  actual  black  boxes  are  composed  of  a  mixture  of  varied 
damping  types. 

TABLE  II  TYPICAL  n  VALUES 


DAMPING 

TYPE 

r\ 

LINEAR 

1.0 

INTERNAL  HYSTERESIS 

0.833 

NON-LINEAR  SPRING 
(INCREASING  STIFFNESS) 

*1.2 

COULOMB  FRICTION 

1.5 

a 

i 


Table  II  shows  that  n  ■  0.833  for  internal  hysteresis.  Consider  the  case 
of  a  solder  joint  on  a  multi-layer  board  (MLB).  There  exists  some  amount 
of  internal  hysteresis  damping  due  to  the  solder  joint  stressing.  However, 
it  is  more  likely  that  the  relationship  between  solder  joint  stress  a  and 
x  of  equation  (4)  is  governed  by  the  friction  damping  at  the  MLB  support 
edges  or  by  some  special  viscoelastic  MLB  damping  treatment.  Therefore,  for 
this  case  in  general  n  i  0.833. 


The  median  stress  cycles  can  be  related  to  test  duration  T  as  follows: 
Nm  “  feff  T  <5) 


fg£f  *  effective  frequency  (Hz) 
rate  of  zero  crossings 


For  a  single-degree-of-freedom  system  [2] 

feff  *  f  M  center  frequency  of  response  spectrum  - 
resonant  frequency  (6 


For  a  two-degree-of-freedom  system  (2DF)  [3] 

r 2 - 2 — p 

t  -  /  a  tr  +  b  fv 


■  1st  mode  resonant  frequency  (Hz) 

■  2nd  mode  resonant  frequency  (Hz) 

■  1st  node  rns  stress  level  (ksi) 

*  2nd  node  rms  stress  level  (ksi) 

■  total  rms  stress  level  (ksi) 
f a  <  feff  "  fb 


W ~  **> 


"  <-  , 


'  *  v,xyn  > rr 


txf'-v- 


0 


T 


(8) 


The  tvo-degree-of-freedom  system  will  later  be  used  in  an  example. 


ACCELERATION  FACTOR  EXPRESSIONS 

Two  separate  acceleration  factor  expressions  will  be  presented.  The  first 
expression  is  for  usually  thought  of  fatigue  case  where  any  existing  initial 
cracks  (i.e.  flaws)  are  not  considered.  The  second  expression  is  the  Fracture 
Mechanics  case  where  initial  cracks,  either  actual  or  postulated,  are 
considered  to  exist.  Both  expressions  are  of  practical  importance.  The 
service  environment  parameters  will  be  denoted  by  the  subscript  1.  The 
accelerated  test  environment  parameters  will  be  denoted  by  the  subscript  2. 


NO  INITIAL  CRACKS 


Reference  [1]  shows  that  the  appropriate  acceleration  factor  expression  is 

"  £oll”,:  ,  ,1/ne 

x,  l-=~)  (g  rms)  (9) 


1  w 


EXAMPLE  1 

Given:  The  service  parameters 

x-,  »  1.0  g  rms;  n  *  0.833 

T^  *  1000  hours 

The  desired  T-,  *  1  hour 


Find:  x-, 

Solution: 

•  • 
x2 


for  a  black  box  witn  copper  wire  as  the  critical  structural  element. 
From  Table  I  S  ■  9.28 


.  i.o  (im)w-n .  z. 


44  g  rms 


1-8 


DUCTILITY  EFFECT 


Equation  (9)  can  be  used  to  show  the  e£fect  of  the  material's  ductility. 
Ductility  is  the  ability  of  the  material  to  be  deformed  without  fracturing. 
Ductile  materials  have  values  of  S  *  9.  Brittle  materials  have  values  of 
8  »  22.  See  Table  1.  Equation  (9)  can  be  plotted  in  a  normalized  fashion. 
Figure  2  shows  a  plot  of  acceleration  for  two  diverse  6  values  with 
n  »  0.833.  It  can  be  seen  that  the  acceleration  factor  is  sensitive  to  8 
values  (l.e.  ductility);  brittle  materials  being  the  most  sensitive. 

DAMPING  LINEARITY  EFFECT 

Figure  3  is  a  plot  of  equation  (9)  for  two  n  values  with  8  -  9.28.  For  a 
time  compression  factor  of  1000  the  acceleration  factor  is  2.44  for 
n  *  0.833  and  1.86  for  n  »  1.2.  The  acceleration  factor  is  sensitive  to 
n  values;  large  n  values  being  the  most  sensitive. 

EQUAL  DAMAGE 

Equal  cumulative  fatigue  damage  exists  at  all  points  along  any  one  line 
(l.e.  curve)  of  figures  2  and  3.  However,  the  damage  is  not  equal  from 
curve  to  curve,  even  at  the  point  where  all  the  curves  intersect. 


o 


o 


FIGURE  2  Acceleration  Factor  for  Variable  B 
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FIGURE  3  Acceleration  Factor  for  Variable  n 


FRACTURE  MECHANICS  EFFECTS 


Initial  cracks  (flaws)  can  exist  in  structural  elements  as  a  result  of 
metallurgical  inclusions,  of  fabrication  or  assembly  procedures,  or  temporary 
overloads.  This  amounts  to  damage  being  cumulated  prior  to  the  service  or 
test  environments.  The  acceleration  factor  expression  is  as  follows: 


(10) 


(ID 


C 


5 


(ksi/g  rms) 


(12) 


a^  ■  initial  crack  length  (inches) 

Y  ■  geometrical  parameter 

AKC  »  material's  fracture  toughness  (ksi  ♦'in  ) 

6  »  slope  parameter  of  material's  crack  growth 

race  curve. 


Equation  (10)  is  a  transcendental  equation  that  is  most  conveniently  solved  by  a 
method  developed  in  reference  [1].  It  cannot  be  normalised  like  equation  (9) 
because  of  the  inherent  non-linearity  of  the  Fracture  Mechanics  process. 
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Table  III  shows  several  typical  0  values.  By  comparing  Table  III  with 
Table  I  It  can  be  seen  that  0*8/ 2.  On  a  log  -  log  plot  the  curve  of 
equation  (9)  Is  (l/n8)  whereas  the  slope  of  equation  (10)  Is  (l/ri6) •  This 
fact  makes  the  acceleration  factor  of  equation  (10)  quite  different  from 
that  obtained  by  equation  (9)  due  to  the  mere  existence  of  Initial  cracks. 
Compare  Figure  4  of  Example  2  below  with  Figures  2  and  3. 


TABLE  III  TYPICAL  0  VALUES 


MATERIAL 

e 

Cr-Mo-V 

4.09 

4340 

4.65 

7073-T6 

4.00 

The  fatigue  life  of  a  structural  element  is  greatly  reduced  by  the  existence 
of  initial  cracks.  However,  the  acceleration  factor  obtained  from  equation 
(10)  is  not  very  sensitive  to  various  values.  It  is  sensitive  to  n  values. 

EXAMPLE  2 

Given:  Critical  element  material:  707S-T6  Aluminum  alloy 
a^.  *  0.007  inches 

x*l  »  1.0  g  rms;  C^  ■  1.0  ksi/g  rms 
T^  ■  1000  hours 

Find:  x*2  as  a  function  of  T  (0.5  hours  _<  T2  £  1000  hours) 

for  n  ■  0.833,  1,  1.2. 


Solution:  Use  equations  (10)  -  (12).  The  results  are  shovn  in  figure  4. 

It  can  be  sees  that  the  acceleration  factor  is  sensitive  to  the 
value  of  n.  For  Tj  ■  1  hour  x^  a  S  g  rms  for  n  ■  0.833. 
x*2  “  2.4  g  rms  in  example  1  which  has  no  initial  crack.  Thus 
Fracture  Mechanics  effects  greatly  Influence  the  acceleration 
factor  value.  As  previously  mentioned  in  the  EQUAL  DAMAGE  section 
the  cumulative  fatigue  damage  is  the  same  along  all  points  of  a 
given  curve  but  not  between  curves. 


MULTI-FACTOR  CONSIDERATIONS 

In  general  all  the  structural  elements  in  a  relatively  complex  electronic 
black  box  will  not  have  the  same  acceleration  factor.  Yet  a  single  factor 
value  must  be  chosen  for  the  accelerated  test.  Such  a  selection  is 
considered  to  be  subjective.  The  value  chosen  will  result  in  a  proper  test 
for  only  one  class  of  structural  elements.  The  other  elements  will  be 
either  under  or  over -tested.  An  average  factor  would  give  average  results. 
The  most  conservative  approach  would  be  to  select  the  largest  acceleration 


factor  value. 


I 
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SIMILITUDE  VIOLATION  COMPENSATION 


Some  similitude  violations  can  be  compensated ' 
be  shown. 


One  such  example  will 


EXAMPLE  3 

Given:  The  similitude  violation  is  due  to  a  difference  in  stress  spectra 

between  service  and  test  environments.  The  stress  vibration  system 
is  2DF.  The  calculated  acceleration  factor  (jL/x^)  is  3  for  the 
desired  time  compression  factor  O^/T^)  using  equation  (9).  The 
spectra  parameter  values  are  given  below: 


PARAMETER 


oa  (ksi) 
(ksi) 


fa  (Hz) 
fb  (Hz) 
oT  (ksi) 
f.f£  <B«) 


ENVIRONMENT 

SERVICE 

TEST 

8 

24 

16 

20 

150 

150 

375 

375 

17.9 

31.2 

342 

264 

Find:  the  appropriate  test  compensation  factors  such  that  oT  *  3  uT 

iTEST  SERVICE 


and  t  ,,  ■  f  .. 

effTEST  eff SERVICE 


Solution:  It  can  be  seen  that  the  resonant  frequencies  are  the  same  at 


both  environments.  However,  o  did  not  increase  from  16  to  48  ksi 

D 


as  desired.  This  would  cause  an  inappropriate  test  damage  state 

•  • 

and  rate,  a-  should  be  3  x  17.9  ■  53.7  ksi.  Therefore,  x„ 
i-rrsT  2 


I 


needs  eo  be  increased  by  an  additional  factor  of  53.7/31.2  » 
1.72.  The  test  duration  needs  to  be  increased  from  its  computed 
compressed  value  by  a  factor  of  342/264  -  1.3. 


EXAMPLE  4 

Given:  The  similitude  violation  is  due  to  a  difference  in  the  amplitude 
distribution  of  the  stress  peaks.  The  motion  of  the  structural 
element  being  stressed  vill  be  snubbed  (i.e.  limited)  if  the  motion 
exceeds  a  specified  displacement.  At  the  service  vibration  level 
motion  limiting  occurs  such  that  the  stress  is  limited  at  5o.  At 
the  accelerated  test  level  the  stress  is  limited  at  3o.  The 
structural  element  being  stressed  is  copper  wire,  xt,  was  determined 
for  the  desired  time  compression  factor  (T0/T^)  using  equation  (9). 

T2  was  computed  to  be  30  minutes. 

Find:  The  modified  time  T2  such  that  equal  damage  is  done  at  both  service 

and  test  levels. 

Solution:  Most  fatigue  damage  during  random  vibration  is  caused  by  stress 
peaks  between  2c  and  5c.  Reference  [1]  shows  the  following  for 
copper  wire: 

1.  Limiting  stresses  at  5c  is  equivalent  to  no  limiting. 

2.  Limiting  stresses  at  3c  will  extend  the  fatigue  life  by  a 
factor  of  1.86. 

3.  5c  peaks  occur  approximately  every  9  minutes  on  the  average. 

5c  peaks  would  occur  during  the  30  minute  test  duration  if  motion  were  not 

limited.  The  compensation  technique  is  to  extend  the  test  time  by  a  factor  of 

1.86.  Thus,  7-,  •  1.86  x  30  ■  56  minutes. 
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